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Introduction 


The  research  project  "Rock  Mass  Persistence"  funded  by  ARO  under  Grant 
No.  DAAG-29-83-K-001  6  started  December  1,  1982  and  ended  November  30,  1985. 
Its  objectives  were  to  first  develop  analytical  models  which  can  describe 
joint  geometries  as  they  occur  in  real  rock  masses.  The  analytical 
descriptions  have  to  be  derivable  from  data  as  obtained  in  usual  field 
exploration.  These  geometric  models  were  then  to  be  used  to  develop  methods 
for  describing  flow  through  jointed  rock  masses  and  for  describing  the 
fracturing,  i.e.,  the  strength  and  def ormability  of  jointed  rock  masses. 

The  objectives  were  achieved.  Basic  knowledge  on  jointed  rock  mass 
description  and  behavior  was  expanded.  Practically  applicable  methods  were 
developed.  The  basis  for  future  fundamental  and  applied  research  was 
created. 

In  addition  to  the  originally  defined  research  objectives  significant 
additional  work  was  performed.  This  consisted  of  the  final  development  of  a 
coupled  flow  deformation  analysis  method  based  on  the  distinct  element 
procedure  and  particularly  of  its  use  in  parametric  studies. 

This  final  technical  report  presents  the  research  results  in  three 
volumes: 

Vol.  I  Rock  Joint  Systems 

Vol.  II  Coupled  Flow  -  Deformation  Analysis 

Vol.  Ill  Fracture  Propagation  in  Jointed  Media 

Each  volume  was  developed  from  a  P'n.D.  thesis.  Given  their  extent 
the  most  important  results  are  summarized  in  this  Executive  Summary.  Each  o: 
the  following  Chapters  2  to  4  is  therefore  devoted  to  a  summary  of 


Volumes  I  to  III  respectively.  This  is  followed  by  general  conclusions  and 
recommendations  for  further  research,  a  listing  of  personnel  who  have  been 
working  on  the  project  and  acknowledgements.  A  listing  of  papers  which  have 
been  or  will  De  produced  Dased  on  this  research  will  conclude  the  Executive 
Summary . 


2.  Rock  Joint  Systems 


2.1  Background 

Geometric  and  mechanical  characterization  of  rock  3omts  is  the  basis  for 
most  of  the  work  of  engineering  geologists,  and  civil  and  mining  engineers 
when  dealing  with  rock  masses.  Also,  such  a  characterization  plays  an 
important  role  in  investigations  on  ;joint  genesis.  Joints  are  usually 
characterized  by  the  parameters  listed  in  Table  2.1;  in  addition  joint 
termination  and  autocorrelation  have  been  considered  in  this  research. 

However,  given  their  pervasive  and  three-dimensional  nature  on  the  one  hand 
and  their  limited  exposure  in  outcrops,  borings  and  tunnels  on  the  other  hand 
makes  complete  descriptions  of  joints  difficult.  Such  an  ideal 
characterization  would  involve  the  specific  description  of  each  joint  in  the 
rock  mass,  exactly  defining  for  instance  its  location,  shape,  size,  planarity, 
aperture  and  shearing  resistance.  This  is  not  possible  for  a  number  of 
reasons:  1)  the  visible  parts  of  joints  are  limited  (for  instance  to  joint 

traces  only),  thus  preventing  a  complete  description;  2)  joints  at  a 
distance  from  the  exposed  rock  surfaces  cannot  be  directly  observed;  3)  direct 
(visual  or  contact  measurements)  and  indirect  (geophysical)  observations  have 
limited  accuracies. 

For  these  reasons  joints  in  a  rock  mass  are  usually  described  as  an 
assemblage  rather  than  as  individual  features.  The  assemblage  has  stochastic 
character  in  that  joint  characteristics  vary  in  space.  Such  variations  may  be 
minute  as  in  the  case  of  the  orientation  of  a  set  of  approximately  parallel 
joints  or  they  may  be  large  if  a  particular  property  has  substantial 
variability.  It  is  important  to  note  that  spatial  variability  can  but  does 


Parameter 


Descnptionn 


Attitude  (orientation)  Strike  and  Dip,  or  Dip  Direction  and  Dip  or  Pole 

Orientation  with  spherical  coordinates. 

Location  Coordinates  of  points  on  iomt  planes  (usually 

spacing  and  location  is  measured). 

Spacing  Distance  between  intersection  points  produced  by 

intersecting  a  line  perpendicular  to  the  joint 
mean  attitude  and  the  loints  of  a  set.  A  set  of 
joints  consists  of  joints  that  are  parallel  or 
subparrallel  to  each  other. 

Size  Extent  of  jointing  generally  expressed  as  length 

of  traces  on  two-dimensional  outcrops  or  as  the 
area  of  jointed  segments  in  a  joint  plane 
(see  also  persistence). 

Persistence  The  proportion  of  actual  jointed  parts  a1  of  a 

joint  plane  A, 


Aperture  Width  of  open  joint  measured  perpendicularly  to 

joint  walls  (only  joints  filled  with  air  or  wate 
are  considered  open,  while  joints  containing 
filler  between  the  wall  rock  are  considered 
closed . 

Joint  Filler  Material  weathered  from  wall  rock,  material 

sedimented  in  joint,  material  precipitated  in 
joint. 

Planarity  Three  dimensional  shape  of  joint  surface; 

generally  either  planar,  or  with  regular  or 
irregular  deviations  from  a  plane. 

Roughness  Deviation  from  the  ideal  planar  surface  on  the 

microscale  usually  witn  wavelength  <  1  cm. 

Waviness  Deviation  from  the  ideal  planar  surface  on  the 

macroscale  usually  with  wavelenctns  >  1  cm. 

Shape  Shape  of  joint  boundaries,  generally  either 

circular,  elliptical,  polvaona!  or  irreauiar. 
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not  nave  to  imply  random  underlying  mechanisms;  spatial  variability  may  gust 
as  well  oe  the  result  of  a  number  of  simultaneous  or  sequential  Dut  very 
specific  mechanisms. 

Two  major  approacnes  nave  emerged  to  describe  the  assemblage  of  joint 
characteristics  in  a  rock  mass:  the  traditional  disaggregate  characterization 
and  the  more  recent  aggregate  characterization .  In  the  former,  each  joint 
characteristic  is  described  separately,  for  instance  through  orientation 
distributions  (pole  diagrams),  spacing  distributions  and  others.  In  the 
latter  the  interdependence  of  joint  characteristics  is  captured  through  the 
formulation  of  joint  system  models.  A  particular  joint  system  model 
represents  a  typical  geometry.  The  individual  characteristics  are  still 
stochastic  but  their  interdependence  is  specified;  for  instance,  in  an 
orthogonal  model  the  mean  orientations  of  three  joint  sets  are  at  right  angles 
to  each  other,  but  some  variability  of  orientation  about  these  means  exists. 

Table  2.2  summarizes  the  rock  joint  system  models.  Each  of  the  five 
models  consists  of  a  particular  combination  of  the  rock  joint  system 
characteristics  m  Table  2.1.  Joint  planarity  is  specified  as  planar  for  all 
joint  system  models,  and  any  jcint  location  or  autocorrelation  process  is 
permissible  within  any  model.  Any  component  specified  as  stochastic  may  also 
be  deterministic. 

The  orthogonal  model  was  developed  by  Irraay*  (1955),  Childs  (1957)  and 
Snow  (1965).  They  applied  the  model  as  did  Schwartz  et  ai.  (1980)  to  hydro- 
logic  problems.  Figure  2.1  shows  a  regular  orthogonal  joint  system  model. 

The  Baecher  model  (Baecner  et  al.  1978)  represents  joints  as  circular  or 

Literature  references  are  listed  in  Appendix  I.  Note,  however,  that  only  a 
very  limited  number  of  references  are  used  in  this  Executive  Summary.  Complete 
literature  references  are  given  in  Volumes  I  to  III. 
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Table  2.2 


JOINT  SYSTEM  MODELS 


MODEL  NAME 

JOINT 

SHAPE 

JOINT 

SIZE 

TERMINATION 

AT  INTERSECT. 

CO¬ 

PLANARITY 

ORIENTATION 

OF  SETS 

1 .  Orthogonal 

Rectangle 

Bounded 

Unbounded 

Unbounded 

no 

yes 

no 

yes 

yes 

Parallel 

Parallel 

Parallel 

2.  Baecher 

Circle 

Ellipse 

Bounded 

no 

no 

Stochastic 

3.  Veneziano 

Polygon 

Bounded 

in  joint 
planes  only 

yes 

Stochastic 

4.  Dershowitz 

Polygon 

Bounded 

yes 

yes 

Stochastic 

5.  Mosaic 

Tesselation 

Polygon 

Bounded 

yes 

yes 

Regular 

Stochastic 

For  all  models,  joints  are  planar,  and  any  location  or  autocorrelation  process 
xs  possible.  Joint  locations  are  usually  stochastic.  Bounding  of  loints 
implies  that  joints  which  are  smaller  than  the  region  under  consideration  can 
De  represented.  This  means,  with  reference  to  the  sketch  in  Table  1,  that  A  is 
tne  region  anc  a^  are  bounded  loints;  if  the  joints  were  unbounded  a^  would  be  > 
A.  Joint  sizes  are  usually  stochastic  either  specified  directly  or  indirectly 
through  stochastic  location  or  orientation. 


elliptical  disks  (Figure  2.2).  The  joint  radii  can  be  regularly  or 
stochastically  distributed  (lognormal  or  exponential  distribution).  Similarly, 
joint  orientations  can  be  regular  (e.g.  exactly  parallel)  or  represented  by  any 
orientation  distribution  (Fisher,  bi-variate  Fisher,  Bingham,  uniform,  or 
bivariate  normal);  also  joint  location  can  be  regular  or  stochastic  (usually  by 
a  Poisson  process) . 

The  Veneziano  model  (Veneziano  1979)  describes  joints  as  polygons.  It  is 
based  on  Poisson  plane  and  Poisson  line  processes  as  shown  in  Figure  2.3. 
Orientation  distributions  are  as  in  the  Baecher  model. 

The  Dershowitz  model  is  also  based  on  Poisson  plane  and  Poisson  line 
process  but  such  that  joint  plane  intersections  produce  joints  (Figure  2.4). 

The  Dershowitz,  Baecher  and  Veneziano  models  have  been  developed  in  this  and 
previous  MIT  research. 

A  variety  of  two-dimensional  and  three-dimensional  Mossaic  Tesselation 
models  (Figs.  2.5  and  2.6)  represent  other  possible  joint  system  models. 

2.2  Joint  System  Modelling  and  Simplified  Expressions 

Joint  System  models  can  be  expressed  by  numerical  procedures  which 
represent  the  two  or  three-dimensional  deterministic  (regular)  and  stochastic 
character  of  the  particular  model.  The  computer  code  JINX  has  been  developed 
for  this  purpose.  It  includes  the  Orthogonal,  Baecher,  Veneziano  and  Dershowitz 
models  with  the  corresponding  deterministic  and  stochastic  characteristics  in 
two  and  three  dimensions.  It  is  thus  the  first  three-dimensional  stochastic 
joint  geometry  code.  JINX  consists  of  6  basic  components  which  can  be  linked 
together  in  any  combination  as  desired.  Each  of  the  components  consists  of  one 
or  more  programs  for  specific  applications.  The  components  are  as  follows: 

•  Common  Variable  Storage:  For  reallocation  of  storage  for  different 


models. 


M 


b)  Poisson  Line 
Process  Formed 
by  Intersections 


c) 

I 


Marking  of 
Polygonal  Joints 


Figure  2,4  Generation  of  Dershowitz  Joi 


Orthogonal  Hexagonal  Triangular 

a.  Regular  (Determinist i c) Tesselat ions 


Poisson  Planes  Polvhedra 

b.  Stochastic  Tesselations 


Figure  2.6  Three  Dimensional  Mosaic  Tesselation 


•  Input  and  Control:  Controls  execution  of  simulation ,  analysis,  anc 
finite  element  preprocessor  packages. 

•  Simulation:  Generation  of  baecner  model,  Veneziano  model  and  Ders'nowitz 
model.  The  usual  point  cnaractens tics  as  listed  in  Table  1  serve  as 
input  parameters. 

•  Analysis:  Calculation  of  joint  intersections  and  sample  statistics 

( RQD,  fracture  frequency)  for  two  and  tnree  dimensional  networks,  of 
point  network  connectivity  measures  and  definition  of  connected  joint 
units  for  two  and  three  dimensional  models.  Calculation  of  block  size 
measures  for  two  and  three  dimensional  networks.  Statistics  and  output 
for  Monte  Carlo  simulations. 

•  Graphics:  Graphical  display  of  two  dimensional  rock  point  networks  and 
projection  of  an  individual  3-D  Dersnowitz  model  joint  plane  on  the 
horizontal  (x-y)  plane.  Stereoplotting  of  generated  point  orientation 
distributions . 

o  Finite  element  flow  programs  and  preprocessors  (see  Section  2.4, 

The  computer  program  package  JINX  is  of  great  significance: 

©  Practically,  it  serves  as  a  basis  for  hydraulic  flow  modelling  (see 
Section  2.4)  in  jointed  media  and  for  fracture  propagation  models 
of  jointed  media. 

•  It  is  an  essential  research  tool  for  the  comparisons  conducted  ir.  this 
research  (see  below)  and  for  future  work  on  joint  genesis. 

Comparisons  between  a  large  number  of  joint  maps  ar.d  photographs  were  made. 
They  showed  that  the  different  joint  system  models  are  necessary  to  cover 
the  variety  of  real  joint  geometries.  While  the  joint  system  models  represent 
a  wide  range  of  real  geometries  they  still  do  not  cover  all  of  them.  Notably, 
they  do  not  cover  geometries  with  non  planar  joints. 


Tne  stocnastic  representations  involve  simulation  with  a  nur.per  of 
realizations  of  tne  particular  model.  It  is  often  aesiraole  for  engineering 
surposes  but  also  in  many  research  applications  to  express  joint  geometry  with 
simpler  measures  ratner  than  to  use  the  computer  simulation.  Extensive  work  in 
this  direction  was  perfomea: 

o  A  numDer  of  so-called  persistence  measures  were  created.  Tney 
involve  expressing  joint  density  as  joint  area  per  volume  or  by 
numDer  of  joints  per  volume. 

•  Several  so  called  connectivity  measures  were  also  developed.  These 
involve  geometric  expressions  which  are  related  to  hydraulic 
conductivity  of  the  jointed  medium  or  which  directly  express  this 
conductivity  (see  also  Section  2.4). 

•  Finally,  several  block  size  measures  were  created  in  which  block  shape 
and  size  are  expressed  m  a  simplified  manner. 

"hese  simplified  measures  have  peen  developed  Doth  for  two  and  three 
dimensions.  Wherever  possible  analytical  expressions  were  derived,  otherwise 
the  measures  were  numerically  ootamed  by  performing  a  number  of  simulations 
witr,  the  complete  moae_s. 

2 . 3  Comparisons 

Extensive  parametric  studies  were  conducted  to  compare  the  different  models 
A  numDer  of  interesting  results  which  were  obtained,  will  De  summarized  here.  By 
running  simulations  with  the  complete  models  and  simultaneously  deriving  the 
simplified  measures  one  can  assess  now  represenati ve  each  of  tnese  measures  is. 
This  makes  it  possible  to  select  tne  more  representative  measures  for  engmeenn 
application,  and  further  research.  From  tnese  investigations  one  can  conclude 

i  c 


that  the  persistence  measure  P2  (P32  in  three  dimension) ,  which  is  the  joint 
trace  length  per  area  (joint  area  per  volume),  the  connectivity  measures  C7 
(which  is  the  equivalent  hydraulic  conductivity)  and  Cg  (equivalent  joint 
network  length,  see  Fig.  2.7)  provide  the  best  simplified  measures. 

Another  interesting  result  is  the  importance  of  joint  connectivity.  The 
Baecher  and  Veneziano  models  with  essentially  independent  joints  and  the 
Dershowitz  model  with  joint  termination  at  joint  intersections  provide  extremes 
for  such  a  comparison. 

The  models  were  compared  in  both  two  and  three  dimensions  using  analytical 
techniques  and  computer  simulations.  Figure  2.8  illustrates  the  relationship 
between  intersection  intensity  C-|  (number  of  joint  intersections  per  unit 
volume)  and  persistence  measure  P32  (area  of  joints  per  unit  volume)  for  disk 
(Baecher)  and  Polygon  (Dershowitz)  joint  system  models.  Clearly,  the 
intersection  intensity  of  the  Dershowitz  model  is  higher  than  that  in  the 
Baecher  model. 

This  difference  can  be  seen  in  terms  of  engineering  parameters  in  Figure 
2.9  which  plots  effective  hydraulic  conductivitv  K  against  joint  size  for 
two-dimensional  models.  Substantial  flows  occur  at  much  lower  levels  of 
intensity  in  the  Dershowitz  model  than  in  the  Baecher  model.  Figure  2.10  shows 
a  similar  result  for  three-dimensional  hydrological  modeling.  In  addition, 
this  figure  shows  that  the  behavior  is  significantly  different  in  two 
and  three  dimensions. 

This  points  out  the  importance  of  both  the  geometry  of  the  joint  system 
model  and  the  dimensionality  of  modeling:  The  use  of  two  rather  than 
three-dimensional  models  may  have  dramatic  effects  on  predictions  of  engineering 


behavior. 


(A)  all  joints  form  the  "joint  system" 

(B)  two  or  more  interconnected  joints  form 
a  "network” 

(C)  joints  with  no  intersections  are  "isolatec 
j  oints" 

CD)  joint  networks  and  isolated  joints  are 
"joint  units" 

Note:  C&  can  be  in  any  direction.  In  this  figure, 

Cgx is  the  x  direction  as  shown. 


igure  2,7  Equivalent  Joint  Network  Length  Cc 
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Figure  2.9  Relation  Between  Effective  Hydraulic  Conductivity 
and  Joint  Size:  Two  Dimensional  Riddels 
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2.4  Hydrologic  Applications  of  Joint  System  Models 

By  combining  tne  joint  geometry  models  whit'n  a  suitable  model  representing 
flow  in  jointed  media  one  can  determine  the  hydraulic  oerformance  of  such 
media. 

The  flow  of  water  through  a  system  of  joints  is  determined  when  the  head 
along  each  individual  joint  is  known  and  such  that: 

a)  Laplace's  equation  is  satisfied  within  each  joint. 

b)  There  are  no  jumps  in  total  head  along  the  intersections  of  two  or  more 
joints  (i.e.,  where  two  or  more  jonts  meet  they  nave  the  same  head). 

c)  Continuity  is  satisfied  along  the  intersections  of  two  or  more  joints 
(i.e.,  the  sum  of  the  influxes  into  two  or  more  intersecting  joints  is 
zero  along  any  segment  of  their  line  of  intersection). 

d)  Tne  conditions  of  no-flow  through  the  free  surface  and  zero  piezometer 

head  (assuming  the  atmospheric  pressure  to  be  zero)  at  the  free  surface 

are  satisfied. 

Such  a  distribution  of  heads  throughout  a  system  of  joints  can  be  found  using  a 
Finite  Element  approach. 

A  computer  program  for  three-dimensional  laminar  flow  through  jointed  rock 
was  developed.  It  uses  Finite  Element  discretizations  of  joints  such  that  the 
Finite  Element  meshes  of  two  or  more  intersecting  joints  share  all  nodes  that  lie 
on  the  line  of  intersection.  The  unknowns  in  this  formulation  are  the  heads  at 
the  nodes.  The  head  at  any  point  inside  an  element  is  obtained  by  interpolation 

from  the  heads  at  the  nodes  of  the  element.  Laplace's  equation  is  satisfied 

inside  each  element  for  any  values  of  the  head  at  the  nodes,  given  an 
appropriate  selection  of  element  type  (e.g.,  four  node  quadrilateral  elements). 
The  heads  at  the  nodes  are  determined  from  the  imposition  of  the  requirement  of 


continuity  at  the  edges  of  the  elements.  If  the  problem  involves  a  free 


surface,  iterations  are  necessary  to  determine  the  position  of  tne  free  surface 
Revising  tne  Finite  Element  mesh  m  each  iteration,  as  is  commonly  done  for  the 
solution  of  two-dimensional  problems  of  free-surface  flow  through  porous  media, 
is  inconvenient  Decause  nere  the  elements  are  tied  to  tne  geometry  of  tne 
joints.  A  scheme  that  does  not  involve  mesh  revisions  is  therefore  used.  The 
Finite  Element  joint  flow  model  and  the  geometric  joint  system  models  were 
combined  in  the  computer  program  package  JINX.  This  was  then  used  to  develop 
some  of  the  connectivity  measures  (Section  2.2  and  2.3)  and  to  develop  a  number 
of  hydrological  applications.  These  applications  can  be  summarized  as  follows: 
e  Scale  effect:  Groundwater  can  flow  in  individual  joints  only  as  far 

as  the  size  of  the  joint.  This  results  in  a  scale  effect  in  which 
hydraulic  conductivity  is  proportional  to  the  intergral  of  the  joint 
size  distribution  over  all  joints  which  are  larger  than  twice  the  scale 
of  interest  (Figure  2.11).  When  the  intersection  of  joints  to  form  join- 
networks^  is  considered,  it  is  the  size  distribution  of  these  networks 
(Cg)  rather  than  the  size  of  individual  joints  which  controls  flow  at 
any  scale  (Figure  2.7).  Figures  2.9  and  2.10  showed  effective  hydraulic 
conductivity  for  networks  of  disk  (Baecher)  and  polygon  (Dershowitz) 
joint  system  models.  For  moth  models,  the  effect  of  scale  is 
pronounced:  As  the  scale  of  the  problem  under  consideration  increases, 

the  effective  hydraulic  conductivity  rapidly  decreases. 

The  result  shown  in  Figure  2.9  and  2.10  has  practical  applications  for 
performance  modeling  of  jointed  rock.  A  major  objection  to 
the  use  of  flow  models  is  the  difficulty  of  measurement  of  effective 
hydraulic  apertures  for  joints.  Tne  shape  of  the  curve  in  Figure  2.9, 
however,  can  be  established  from  statistics  of  icir.t  orientation,  trace 


Probab 


length,  and  intensity,  which  can  readily  oe  ontaine;  from  oorenoies  and 
napping  in  outcrops  and  adits.  Tne  actual  values  on  tne  ordinate 
(effective  hydraulic  conductivity)  can  tnen  oe  oDtained  by  oerforminc 
nydroloaical  tests  at  two  or  more  scales  and  tnus  normalizing  tne  curve. 
This  allows  one  to  make  croiections  about  tne  effective  nyaraulic 
conductivity  at  much  larger  scales. 

Anisotropy  Effect:  Just  as  the  distrmution  of  trace  lengths  of 

individual  joints  controls  the  scale  at  wnicn  flow  can  occur  in 
individual  joints,  the  distribution  of  joint  orientation  determines  the 
anisotropy  of  flow  m  individual  joints  (Figure  2.) 2a).  Similarly,  for 
networks  of  interconnected  joints,  the  anisotropy  of  joint  system 
geometry  determines  the  nydrological  anisotropy  (Figure  2.12b). 

Figure  2.13  shows  the  first  quadrant  of  permeability  ellipses  at  scales 
of  2.5,  10  and  20  times  the  mean  joint  length  for  the  Ders'nowitz  joint 
system  model  with  individual  joint  orientations  distributed  according  to 
a  Fisher  distribution  witn  very  low  dispersion.  As  a  result,  at  a  small 
scale,  corresponding  to  the  influence  of  individual  joints,  the  effective 
hydraulic  conductivity  is  highly  anisotropic.  At  larger  scales,  however, 
joint  networks  control  flow  and  the  permeability  ellipse  indicates  much 
less  anisotropy. 

This  result  has  significant  implications  for  tne  assessment  of  the 
hydrological  anisotropy  of  jointed  rock.  The  hydrological  anisotropy  of 
joint  systems  is  scale  dependent  and  depends  strongly  upon  the  formation 
of  joint  networks  of  interconnected  joints.  In  situ  tests  used  to 
evaluate  anisotropy  can  only  be  used  to  predict  anisotropy  at  the  scale 
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scale  of  individual  joints,  anisotropy  of  flow 
directly  related  to  anisotropy  of  joint  orientations 


scale  of  joint  networks,  joint  orientation 
sotropy  may  have  less  effect  on  flow  anisotropy 


Figure  2,12  Scale  Effects  on  Joint  System  Anisotropy 
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at  whicn  they  were  run.  In  order  to  extrapolate  these  results  to  larger 
scales,  joint  system  models  must  De  developed  ana  evaluated  as  in  Figure 
2.13. 

Aperture  Effects:  Many  researchers  (e.g.  Long  et  al.f  1983,  Snow, 

1968)  have  noted  that  since  effective  hydraulic  conductivity  is 
proportional  to  joint  aperture  cuDed  and  the  distribution  of  aperture  is 
generally  assumed  to  be  lognormal,  aperture  uncertainty  and  variability 
tends  to  dominate  all  other  effects  in  determining  joint  system 
hydrological  properties. 

While  this  is  correct  for  the  case  where  all  flow  occurs  in  individual 
joints,  it  is  not  necessarily  so  for  flow  through  networks  of 
interconnectec  joints.  At  large  scales,  where  flow  occurs  in  series 
throuah  interconnected  joints  with  independent  apertures,  the 
variability  of  aperture  has  much  less  influence  on  hydrological 
properties  than  would  be  assumed  based  upon  flow  in  sinale  joints.  This 
result  is  illustrated  by  the  distribution  of  effective  hydraulic 
conductivities  for  flow  through  joint  systems  of  n  joints  with  a 


lognormal  aperture  distribution  (Figure  2.14). 


Effective  Hydraulic  Conductivity,  K  (cm/sec) 


3.  Couoied  Flow-Deformation  Analysis  for  Jointed  Media 


3 . 1  background 

In  many  engineering  problems  in  ]omten  rock,  deformation  and,  eventually, 
instability  (failure)  involve  interaction  Detween  cleft  water  pressure  in  tne 
joints  and  deformation  of  tne  joints.  Closure  or  opening  of  the  joints 
perpendicuarly  to  their  surface  influences  water  flow  and  pressure  in  the 
joints.  Similarly,  increases  and  decreases  of  the  water  pressure  in  joints 
affect  joint  apertures.  This  interaction  Decones  more  complex  if  joints 
displace  along  their  surface  (in  shear);  the  usual  associated  asperity 
interaction  and  dilatancy  effects  can  affect  joint  aperture.  Also,  joint 
fillers  add  to  the  complexity  of  the  problem. 

A  method  of  analysis  which  can  treat  the  interaction  of  joint  deformation, 
cleft  water  pressure,  and  flow  is  needed  when  designing  structures  in  and  on 
rocks  and  also  for  gaining  a  more  basic  understanding  of  jointed  rock  mass 
behavior  from  a  geologic  point  of  view.  The  method  should  also  be  useful  when 
interpreting  monitoring  instrumentation.  For  instance  water  pressure  and 
deformation  in  rock  masses  can  De  measured.  With  an  appropriate  analysis  method 
inferences  on  the  state  of  tne  rock  mass,  particularly  in  relation  to 
instability,  can  be  drawn  from  the  data.  This  will  make  it  possible  to  integrate 
monitoring  into  design  which  is  required  in  the  application  of  advanced 
"observational"  or  "adaptable"  design  methods;  in  such  approaches  the  design  is 
optimally  adapted  to  encountered  conditions,  eliminating  unnecessary- 
conservatism  . 

A  variety  of  continuum  and  aiscontinuum  methods  exist  for  rock  mass 
deformation  and  failure,  and  flow.  In  an  initial  phase  of  the  research  it  was 
decided,  based  on  an  extensive  review  of  available  approaches,  that  ciscontinuum 
methods  had  to  De  used  for  both  the  deformation  anc  the  flow  analysis. 
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Lor.unuur.  approac.nes  are  su.ie"  tc  treatment  cr  ..arut  relative  cisr.  ..acemer.t 

as  tnev  occur  lr.  loints.  Equi  valent  porous  ne.'.;c  ar:  roacne-  are  of  ter. 
unreoresentati ve  wnen  analyzing  jointed  media  a?  was  exriames  n.  Cnanter  1  cf 
this  Executive  Summary.  Consequently  trie  so  called  ncic  r iocr.  or  distinct 
element  method  (Cundall  1971,  1974)  was  cnoser.  tc  monel  aeformatior.  and  it  was 

decided  to  couple  it  with  a  discontinuous  flow  moael. 

3 . 2  Coupled  Flow  Deformation  Analysis 
2.2.1  Deformation  Analysis 

In  the  Rigid  Block  Method,  tne  rock  mass  is  assumed  tc  consist  of  separate 
blocks  touching  eacn  other  at  a  discrete  number  of  contacts  where  elastic, 
frictional  and  viscous  forces  are  applied. 

Two  types  of  initial,  i.e.  before  any  displacement,  bloc):  configurations 
are  possible  in  the  computer  code  developed  here,  as  shown  in  Fig.  3.1.  Tne 
configuration  in  Fig.  2.1a  simulates  two  fully  persistent  joint  sets  of 
arbitrary  (but  constant  for  each  joint)  orientations  and  spacings.  The 
configuration  in  Fig.  3.1b  simulates  two  joint  sets  of  arbitrary  (but  constant 
for  each  joint  set)  orientations  and  spacings,  one  joint  set  being  fully 
persistent  and  the  other  having  a  persistence  of  fifty’  percent.  From  these 
oasic  assemblies  any  shape  of  a  real  jointed  rock  mass,  e.g.,  a  slope  or  the 
rock  mass  around  a  tunnel  can  be  created.  Larger  or  irregularly  shaped  blocks 
are  also  possible,  provided  that  they  can  be  formed  by  combining  elementary 
blocks  configured  as  shown  in  Fig.  3.1. 

In  the  distinct  element  method  the  blocks  are  considered  rigid,  out  tne 
joints  are  considered  deformable,  i.e.  all  deformations  are  represented  as 
joint  deformations.  This  is  a  reasonable  assumption  for  most  rock  mass 
applications.  Deformations  of  the  rock  may  occur  by  relative  displacement  of 
the  rock  clocks  as  snown  m  Fic.  3.2.  Reactions  (forces)  at  tne  contacts  are 


Block  configurations  t 
CAN  EE  COMBINED  TO  FOR 
SHAPED  BLOCKS. 


obtained  from  the  penetration  or  separation  at  contacts  through  the  application 
of  appropriate  deformation  laws.  These  can  include  time  dependent  behavior  or 
the  interaction  with  fillers.  Deformation  analysis  by  dynamic  relaxation 
proceeds  in  time  steps,  evaluating  at  a  series  of  successive  instants  the  forces 
acting  on  each  block  (its  own  weight,  the  reactions  of  neighboring  blocks,  and 
water  forces).  These  forces  determine  the  acceleration  of  each  block,  which  is 
then  integrated  to  find  the  displacement  increment  during  the  time  step  under 
consideration.  This  procedure  is  essentially  that  of  Cundall  (1974)  with  the 
modifications  described  in  3.3.2  and  3.3.3  below. 

3.3.2  Flow  Model 

Flow  is  modeled  as  taking  place  within  a  network  of  conduits  corresponding 
to  block  sides  (which  can  be  thought  as  pipes)  meeting  at  a  number  of  'nodes' 
corresponding  to  block  corners.  In  modeling  flow  the  changes  in  the  joint 
apertures  that  are  produced  by  the  deformation  model  are  taken  into  account. 
However,  except  for  changes  in  apertures,  the  flow  algorithm  assumes  that  the 
geometry  does  not  change  during  deformation.  This  is  equivalent  to  assuming 
that  flow  takes  place  in  a  network  of  conduits  of  varying  cross-sections  but  of 
constant  lengths  and  positions  in  space.  Therefore,  the  geometry  of  this  flow 
network  depends  on  the  initial  position  of  the  blocks,  and  subsequent  changes 
affect  only  apertures.  This  apparent  disparity  between  the  flow  and  the 
deformation  algorithm  is  justified  by  the  fact  that  the  flow  variables  (e.g., 
heads,  pressures)  are  much  more  sensitive  to  aperture  changes  than  to  changes 
in  the  locations  of  nodes. 

The  solution  of  the  flow  problem  is  standard.  Laminar  flow  is  assumed; 
flow  in  joints  with  fillers  described  by  Darcy's  law  and  for  open  joints  by  the 
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(where  q  =  flow  per  unit  width,  a  ■=  aperture,  g  =  acceleration  of  gravity, 
v  =  kinematic  viscosity,  i  =  head  loss  per  unit  length).  Flow  in  joints  is 
obtained  from  heads  at  the  nodes  and  by  satisfying  continuity  conditions  at  the 
nodes.  To  make  the  problem  determinate,  heads  of  some  nodes  need  to  be  given  as 
boundary'  conditions. 

This  procedure  is  complicated  if  a  free  surface  is  involved.  In  the  model 
solution  procedure,  an  assumption  is  made  as  to  which  nodes  are  included  in  the 
flow  region.  Given  the  heads  at  some  of  then,  it  is  possible  to  solve  the  set 
of  equations  obtained  by  applying  the  condition  of  continuity  at  nodes  where 
head  is  unknown  and  thus  to  determine  these  unknown  heads.  Based  on  these  heads 
it  is  possible  to  check,  and  if  necessary,  update  the  assumption  about  which 
nodes  participate  in  the  flow.  After  making  these  corrections,  heads  at  nodes 
included  in  the  flow  are  computed  again.  The  same  procedure  is  repeated  until 
the  assumed  position  of  the  nodes  relative  to  flow  does  not  change. 

While  the  flow  problem  is  solved  iteratively  in  this  manner,  block 
positions  remain  unchanged.  Figure  3.3  shows  how  a  joint's  aperture  used  for 
flow  computations  is  related  tc  the  joint  deformation  as  described  by  the  Block 
Method.  Clearly,  raw  deformation  data  from  the  Rigid  Block  Method  are 
inappropriate  for  studying  flow,  since  the  Rigid  Block  Method  accepts 
overlapping  of  a  joint's  faces.  Therefore,  "rigid  block"  apertures,  which  can 
be  positive  or  negative,  are  converted  to  "real"  apertures,  which  are  always 
positive  and  can  be  used  in  modeling  flow  in  jointed  rock. 

3.2.3  Coupling 

Coupling  is  handled  through  an  iteration  of: 

1.  Obtain  joint  apertures  from  deformation  analysis. 
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Figure  3.3  Relationship  between  Rigid  Block  deformations 

AND  FLOW  APERTURES 


2.  Determine  flow  and  pressure  in  joints  from  flow  model. 

3.  Modify  forces  acting  on  blocks  based  on  Step  2. 

4.  Step  1. 

3.2.4  Computer  Code 

The  Computer  Code  RGDBLK  was  developed  for  the  coupled  flow-deformation 
anlaysis. 

Emphasis  was  placed  on  making  the  computer  program  user-friendly;  it 
incorporates  routines  for  automatically  setting  up  the  geometry  of  a  problem 
(e.g.  a  rock  slope)  and  its  boundary  conditions  for  flow  with  a  minimum  of 
information  from  the  user  (e.g.  joint  orientation  and  spacing,  position  of  the 
water  table,  etc.).  The  program  is  interactive,  making  use  of  a  video  screen, 
and  lets  the  user  intervene  in  the  solution  process.  For  example:  one  can 
model  the  initial  state  of  stress  under  no-flow  (hydrostatic)  conditions 
prevailing  before  a  rock  slope  is  excavated,  followed  by  modeling  the 
"excavation"  of  the  slope  by  removing  some  of  the  blocks,  and  finally,  modeling 
the  coupling  of  flow  and  deformation  as  the  rock  slope  deforms.  The  user  can 
control  the  iterative  solution  procedure  through  several  parameters  such  as  the 
size  and  number  of  time  steps  and  the  degree  of  coupling  between  flow  and 
deformation  (that  is,  the  frequency  with  which  updated  information  about 
deformation  should  be  used  to  update  flow  forces  or  vice  versa) .  The  output  of 
the  program  is  in  the  form  of  plots  of  block  displacements  at  various  instants, 
of  information  about  the  water-head  field  at  various  instants,  and  of 
information  about  forces  and  displacements  at  particular  blocks. 

Due  to  its  considerable  flexibility,  this  computer  program  can  be  used 


either  for  modeling  particular  situations  occuring  in  practice  and  thus  serving 


as  a  design  tool,  or  no  simulate  the  Dehavior  of  jointed  rock  masses  in  a  more 
general  manner,  thus  serving  as  a  research  tool. 

3 . 3  Case  Studies 

The  flow  deformation  model  is  applied  to  three  typical  cases  of  rock  mass 
behavior,  a  rock  slope,  a  well  in  a  jointed  mass  and  a  dam  and  reservoir  on  a 
jointed  mass.  In  all  tnree  cases  extensive  parametic  studies  were  conducted. 

In  the  Sections  below  only  one  typical  example  will  De  presented,  followed  by  a 
summary  of  the  parametric  study  results,  the  details  of  which  are  presented  in 
Volume  IX . 

3.3.1  Slope  in  a  Jointec  Rock  Mass 

The  geometry  of  the  problem  is  shown  in  Fig.  3.4  After  excavation,  the  water 
levels  are  assumed  to  be  m  the  positions  marked  ^A  and  B,  the  position  B 
representing  an  open  crack.  (It  should  be  noted  that  simultaneous  excavation  and 
lowering  of  the  water  table  can  also  be  investigated.)  Tne  phreatic  surface  and 
water  pressure  distribution  between  these  two  points  was  then  determined  using  the 
coupled  model  and  other  approaches.  The  coupled  model  provides  displacements  ir. 
addition  to  flow  results.  In  addition  to  calculating  the  displacements  and  water 
pressure  distributions,  the  latter  were  used  for  a  comparison  study.  "Standard", 
triangular,  trapezoidal  or  rectangular  water  pressure  distributions  are  assumed  m 
the  critical  joint.  The  stability  of  rock  mass  above  this  joint  was  determined 
using  limit  equilibrium  analysis  with  these  pressure  distributions  and  those 
obtained  from  the  coupled  model. 

Three  different  cases  have  been  studied  with  the  coupled  model.  The  first 
two  allow  one  to  compare  the  model  with  other  models.  The  third  makes  full  use  of 
the  coupled  model's  capabilities: 
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1.  Assume  joints  of  identical  aperture  throughout  the  slope,  and  model  flow. 
(This  is  equivalent  to  a  pipe  network  approach). 

2.  Assume  deformation  due  to  excavation  takes  place  and  different  joint 
apertures  result  (smaller  in  zones  of  higher  stresses) .  Flow  is  then 
modeled  through  these  joints  but  no  coupling  is  introduced.  This,  in 
essence,  represents  a  study  that  one  could  do  using  two  separate  models, 
first  one  to  determine  displacements  and  then  a  discrete  flow  model  to 
determine  water  pressures  in  the  deformed  model. 

3.  Assume  coupling  of  apertures  and  flow. 

Figure  3.5  shows  the  deformed  mass  after  excavation,  but  without  flow.  The 
joint  overlaps  indicate  that  aperture  changes  will  have  a  significant  effect  on 
flow  patterns  and  thus  pressures.  Figure  3.6  presents  water  pressure  distributions 
in  the  "critical  joint"  that  were  obtained  in  the  three  above  mentioned  cases 
studied  with  the  coupled  flow  model. 

These  pressure  distributions  should  be  compared  with  the  "usual"  assumptions 
in  Fig.  3.4.  (The  critical  -joint  is  that  used  in  the  standard  limit  equilibrium 
analysis. )  The  pressures  with  the  equal  size  pipe  network  in  case  1  are  greater 
than  with  the  usual  assumptions  (Fig  3.4)  but  smaller  than  in  deformed  joints  in 
case  2  (different  size  pipe  network).  This  pressure  increase  was  to  be  expected 
given  that  joints  close  under  increased  stresses.  When  coupling  is  introduced  in 
case  3  the  water  pressure  tends  to  open  the  joints  and  thus  causes  a  pressure 
decrease  compared  to  case  2.  These  results  are  specific  to  the  case  studied  and 
should  not  be  generalized. 

The  effect  of  these  differences  in  pressure  distribution  on  slope  stability 
was  quite  significant.  The  required  friction  angles  for  the  slope  in  Fig.  3.4  to 


be  stable  are: 


1^0  10  20 

Water  Pressures 

( kN/m2 ) 


Key : 

1.  Joints  with  equal  apertures 

2.  Joints  with  apertures  changed  by  deformation 

3.  Coupled  case 


Figure  3.6  Water  Pressure  Distributions  in  "Critical  Joint" 
for  Cases  Referred  to  in  the  Text 


Triangular  Pressure  Distribution:  34° 

Trapezoidal  Pressure  Distributor  37° 

Distribution  in  equal  size  pipe  network:  39° 

Distribution  in  different  size  pipe  network:  4t° 

Distribution  in  coupled  case:  40° 

These  results  show  that  the  usual  assimptions  may  be  unconservative  in  some  cases 
This  is  especially  true  for  local  instability  where  individial  blocks  become 
unstable. 

A  series  of  parametric  studies  was  also  carried  out  in  which  the  friction 
angle  and  the  joint  parameters  relevant  to  conductivity  of  water  were  varied.  By 
considering  different  friction  angles,  the  appropriateness  of  the  method  for 
predicting  failure  was  established;  moreover,  patterns  of  numerical  behavior 
corresponding  to  stable  or  failing  cases  were  identified.  By  considering 
different  values  for  the  parameters  affecting  conductivity  it  was  possible  to 
study  situations  where  coupling  is  more  or  less  important.  As  might  have  been 
expected,  the  less  uniform  the  apertures  of  the  joint  network  were,  the  more 
pronounced  was  the  effect  of  coupling.  Although  in  real  cases  there  may  be 
considerable  uncertainty  regarding  the  conductivity  properties  of  the  joints,  it 
is  usually  possible  to  estimate  an  upper  bound  of  coupling  influence  for 
particular  cases. 

3.3.2  Well  in  a  Jointed  Rock  Mass 

Injection  or  withdrawal  of  water  in  a  jointed  rock  mass  causes  displacements 
of  blocks,  increasing  and  decreasing  joint  apertures  and  water  pressures.  This  is 
again  a  case  where  coupling  of  flew  and  deformation  seems  to  be  necessary.  Fig. 
3.7  shews  the  displacements  of  blocks  caused  by  an  injection  in  the  center.  Agai 
the  comparison  of  pressures  along  the  center  joint  reveals  important  information. 
Fig.  3.8  shews  pseudo-equipotentials  (also  lines  of  constant  pressure  in  this 
case)  under  the  assumption  of  uniform  apertures  (a),  and  when  changes  in 


Pseudo-Equi potentials  around  Well  in  Network  of  Joints  of 
Equal  Apertures 


apertures  due  to  coupling  are  considered  (b).  Corresponding  flow  rates  are  0.95  x 
10-5  mVsec  and  4.53  x  10-5  m3/sec.  In  this  example,  taking  into  account  joint 
deformations  has  a  pronounced  effect  and  a  coupled  analysis  gives  different 
results  from  an  uncoupled  one. 

As  in  the  case  of  the  slope,  parametric  studies  were  carried  out  varying  the 
parameters  affecting  joint  conductivity.  More  pronounced  coupling  effects 
correspond  to  parameter  values  associated  with  less  uniform  apertures  in  the  joint 
network.  It  was  also  attempted  to  simulate  observed  behavior  of  well  tests.  In 
particular,  the  coupling  of  the  response  of  the  head  at  the  well  to  the 
application  of  flow  rates  that  are  step-wise  functions  of  time  were  simulated. 

This  was  accomplished  with  good  results. 

3.3.3  Dam  and  Reservoir  on  Jointed  Rock  Mass 

Flew  of  water  from  the  reservoir  through  jointed  rock  under  the  dam  (Fig. 

3.9)  contributes  to  water  loss  and  has  therefore  negative  economic  consequences. 
Water  pressure  distribution  under  a  dam  has  significant  effects  on  dam  stability. 
Drainage  and  grout  curtains  are  used  to  control  pressure  buildup  and  flow. 

While  a  number  of  seepage  flow  and  water  pressure  measurements  under  dams 
have  been  conducted,  there  is  a  lack  of  analytical  approaches  for  flow  and 
pressure  prediction.  This  is  true  for  jointed  rock  (see  Chapter  2  of  this 
Executive  Summary),  particularly  for  cases  where  water  pressure  and  joint 
deformation  strongly  interact.  This  situation  has  become  acute  through  the 
recognition  that  in  many  concrete  dams  the  joints  open  up  on  the  upstream  side 
leading  to  pressure  increases  both  at  depth  and  possibly  under  the  dam  itself. 

This  opening  of  joints  is  caused  by  t'  dewnstream  displacement  of  the  dam.  This 
same  displacement  can  also  close  joints  dewnstream  of  the  dam  and  effectively  lead 
to  high  water  pressures  in  the  downstream  ground  with  negative  consequences  on 
stability. 
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Figure  3.9  Dam  Geometry 
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The  dam  parametric  study  was  conducted  to  shed  light  on  these  issues.  The 
following  variables  were  studied: 

•  Existence  (or  non  existence)  of  grout  curtain  or  drainage  curtain. 

e  Effect  of  length  of  these  curtains. 

o  Joint  filler  existence  and  characteristics. 

o  Joint  aperture, 
and  as  before, 

•  Uniform  pipe  network  flw,  deformed  pipe  network  flow,  coupled 
flow-deforma  tion. 

In  Fig.  3.10  and  3.11  a  case  without  drains  or  grout  curtain  is  shown  for  the 
uniform  apertures  (pipe  network)  and  coupled  conditions  respectively.  Clearly 
coupling  effectively  opens  the  upstream  joint  and  this  leads  to  full  reservoir 
pressure  under  the  upstream  comer  of  the  dam.  The  displacements  (Fig  3.10,  3.11) 
reveal  that  such  a  joint  opening  does  indeed  occur,  it  is  also  apparent  in  the 
"uniform  pipe  network"  case  but  is  much  more  pronounced  when  coupling  is 
introduced. 

In  Fig.  3.12  the  no  drain  -  no  grout  situation  is  compared  to  a  grout  curtain 
case  (for  coupled  conditions).  A  grout  curtain  pushes  the  pressures  to  higher 
values  at  depth  where  a  greater  pressure  gradient  occurs.  It  is  basically 
advantageous  to  have  larger  gradients  (and  higher  absolute  values)  at  greater 
depth.  Also  the  pressure  under  the  toe  (tail  end)  of  the  dam  is  less  than  in  the 
no  grout  -no  drain  case.  With  a  drainage  curtain  (not  shown)  the  expected 
complete  pressure  drop  (to  the  tailwater  level)  occurs  at  the  upstream  comer  of 
the  dam,  where  the  drainage  curtain  is  in  contact  with  the  dam. 

While  the  upstream  joint  opening  has  important  effects,  the  closure  of 
dcwnstream  joints  is  such  as  to  not  greatly  affect  the  downstream  pressure 
distribution  and  magnitudes. 
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The  parametric  studies  consisted  of  varying  the  parameters  listed  in  Fig. 
3.13.  Coupling  proved  to  be  significant  in  the  cases  without  drains  (drains 
essentially  isolate  the  dam  foundation  from  the  flew  region) .  The  effect  of 
coupling  was  studied  in  terms  of  the  water  pressure  diagrams  under  the  dam,  the 
factor  of  safety  of  the  dam  as  it  can  be  derived  from  such  diagrams,  and  the  rate 
of  flow  of  reservoir  water  underneath  the  dam.  Ignoring  coupling  was  seen  to  lead 
to  an  underestimation  of  water-pressure  under  the  dam.  In  comparison  to  the 
equal-aperture  assumption,  this  underestimation  reached  approximately  25%,  leading 
If  to  even  more  significant  overestimation  of  the  factor  of  safety.  Estimates  of  the 

losses  of  reservoir  water  were  seen  to  depend  significantly  on  whether  coupling  is 
recognized  and  also  on  the  hard  to  measure  conductivity  properties  of  the  joints. 

3 . 4  Wedge  or  Block  Stability"  Analysis  with  Rotation  and  Time  Effects 
2.4.1  Basic  Considerations 

Failure  of  rock  slopes  often  involves  study  of  wedges  or  blocks  as  shown  in 
Fig.  3.14.  In  the  examples  shown  in  Fig.  3.14  the  wedge  can  either  slide  on  one 
or  two  planes.  Cases  with  blocks  constrained  by -more  than  -two  planes  are  also 
possible.  In  the  one  plane  case  the  motion  of  the  wedge  is  constrained  only  by 
the  requirement  that  the  wedge  stays  on  the  plane  of  sliding.  Such  two- 
dimensional  motion  will,  in  general,  also  involve  rotation  of  the  wedge.  In  the 
two  (or  more)  plane  case  the  wedge  can  move  only  in  the  direction  of  the  inter¬ 
section  of  the  two  planes  on  which  it  slides  and  the  motion  is  one-dimensional. 

In  addition  to  modelling  of  unidirectionl  or  rotational  sliding  it  is  also 
desirable  to  include  time  effects.  The  failure  (instability)  of  rock  slopes 
usually  evolves  over  time  where  "time”  can  be  very  short  (seconds)  or  extend  over 
many  years.  Ideally ,  one  would  like  to  model  wedge/block  instability  in  a  way 


that  allows  one  to  consider  the  entire  range  of  time  dependent  effects. 


JOINTS:  With  filler  -  without  filler 

APERTURES:  small  -  medium  -  hiah 


GROUT  CURTAIN:  absent  -  1  block  deep  -  2  blocks  deep 

3  blocks  deep 

DRAINS:  absent  -  1  block  deep  -  2  blocks  deep  - 

3  blocks  deep 

UPSTREAM  WATER  LEVELS:  low  -  high 

K^:  low  -  high 


igure  3.13  Parameters  Used  in  the  Dam  Case  Studies 


Failure  of  a  wedge  by  sliding  on  two  planes 


Figure  3.14  Wedge  Shapes 


A  rigid  block  type  analysis  and  computer  code  nas  been  developed  to  model 
wedge /block  rotation  with  time  dependence.  The  assumption  of  a  rigid  block  for 
tne  intact  rock  wedge  (block)  Dounaec  dv  gomts  is  a  good  approximation  as 
discussed  in  Section  2.1.  Tne  rigid  block  procedure  is  essentially  tne  same  as 
discussed  in  Section  2.2  witn  the  following  differences: 

•  The  wedge/block  is  a  single  block  Dut  the  goint  surfaces  are  subdivided 
into  elements  (this  will  be  shown  in  2.4.2  Implementation). 

•  A  model  for  time  dependent  goint  benavior  was  developed  (Section  2.4.2) 

2.4.2  Joint  Behavior  Model 

A  model  for  goint  behavior  in  a  strain-controlled  loading  is  postulated  in 
which  shear  strain  increases  monotonicaliy  as  a  function  of  normal  stress,  shear 
displacement  and  rate  of  shear  displacement.  Shear  stress  is  divided  into  two 
additive  components.  Both  compondents  depend  on  the  normal  stress,  one  being  also 
a  function  of  displacement  and  tne  other  of  displacement  and  displacement  rate. 

Tne  second  term  disappears  when  the  displacement  rate  is  small,  but  becomes 
dominant  when  the  displacement  rate  is  large.  Shear  stress  on  the  failure  surface 
is  then  expressed  as 

T  =  ts  (o  ,  6  )  -f  id(c,c,6)  ,  (2.1  ) 

wnere  t  is  shear  stress;  ts  its  'spring'  component,  depending  on  normal  stress 
c  and  shear  displacement  6  ;  its  'dasnpot'  component  depending  on  c,  6 

and  on  shear-displacement  rate,  6  . 

The  'sprinc'  component  of  shear  stress,  (i.e.,  the  rate-independent 
component),  depends  on  shear  displacement  in  a  form  like  that  cf  Figure  2.15a. 
Figure  2.15d  shows  the  assumed  constant-peak-dispiacement  type  cf  dependence  of 
shear  stress  oehavior  on  normal  stress,  while  Figure  2.15c  presents  an  alternative 
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model  displaying  constant  stiffness  behavior .  both  types  of  Denavior  are 
supported  by  examples  of  direct  shear  test  results  in  the  literature. 

Tne  displacements  and  are  parameters  of  the  model.  When  the 

displacement  is  less  than  6j  ,  asperities  of  the  joint  surfaces  deform  without 
breaking  or  sliding.  Tneir  movement  is  constrained  Dy  aspentites  with  which  they 
interlock  and  by  friction  forces  between  the  joint  surfaces.  In  the  region 
0  <  6  <  6j  the  shear  'modulus'  (relating  stress  to  displacement)  of  tne  joint  is 
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tan?  I 


(3.2) 


that  is,  it  increases  with  normal  stress. 

As  long  as  6  <  tne  displacement  is  entirely  due  to  the  deformation  of 

the  asperities.  When  a  displacement  of  6t  is  reached,  the  asperities  are  on  the 
verge  of  creaking  or  sliding.  The  shear  stress  at  this  point  is 


:s  =  c  ~  c  tang- 


13. 3) 


~.e.,  it  is  equal  to  a  component  proportional  to  normal  stress  plus  a  cohesion 
intercept. 

As  tne  displacement  increases  beyond  6-  ,  asperities  start  to  break  or 
slide.  Breaking  is  associated  with  a  drop  in  shear  stress,  anc  the  higher  the 
normal  stress,  the  more  extensive  is  the  breaking  of  asperities  and  the  iaraer 
tne  drop  in  shear  stress.  Under  zero  normal  stress,  sliding  is  predominant  and 
there  is  almost  no  breaking  of  asperities,  hence  no  drop  in  snear  stress.  In  tne 
region  6l  <  6  <  6j-  ,  shear  stress  and  displacement  are  related  by  the  'modulus' 
tan? —  tan?T 


v  =  c 
II  c _ -  c . 


<  0 


(3.4) 


by  tne  time  tne  displacement  reaches  6--  ,  the  breaking  of  asperities  has  ceasec 


as  all  asperities  constraininc  relative  displacement  of  the  two  sides  of  tne  icint 


nave  oeen  eliminated.  For  displacements  larger  tnan  on  ,  tne  shear  stress  is 
given  by 

-  s  =  c  *  c  tampxi  •  (2.5) 

The  'aashpot'  component  of  snear  stress,  Tc  <  is  described  as  a  function  of 
Oj  and  on  in  Figure  3.16.'  It  is  proportional  to  the  normal  stress.  In  the 
range  0<  6  <  oj  ,  Tc2  depends  on  the  rate  of  displacement  6  but  not  on 
displacement  and  is  equal  to  zero  when  c  >  oji  •  This  behavior  is  not  based  on 
experimental  results,  as  creep  along  discontinuities  has  received  very’  little 
attention  until  now.  However  this  assisted  behavior  seems  plausible  if  one 
associates  viscous  behavior  with  deformation  of  intact  rock  that  forms  the 
asperities.  Viscous  effects  should  decrease  in  importance  as  the  relative 
displacement  of  the  two  faces  of  the  goint  becomes  larger,  as  the  displacements 
increase  sliding  over  asperities  or  breaking  of  asperities  rather  than  viscous 
aeformation  of  the  asperities  will  occur. 

2.4.3  Implementation  of  the  Model 

The  computer  program  WEDGE  has  been  prepared  to  apply  the  procedures  outlined 
above.  It  can  accomodate  any  wedge  geometry  and  any  loading  history.  In  its 
present  form  it  computes  normal  stresses  based  on  the  assumption  that  the  wedge  is 
rigid  and  that  the  underlying  rock  has  a  Winkler- type  behavior,  i.e. ,  the  normal 
stress  at  any  point  on  the  failure  surface  is  proportional  to  the  normal 
displacement  at  that  point,  given  that  the  coefficient  of  proportionality  is  the 
same  for  all  points  on  the  failue  surface.  However,  if  desired,  values  for  the 
normal  stresses  as  computed  by  a  Finite  Element  program  can  be  used  as  input.  The 
joint  behavior  model  is  that  cf  Section  3.4.2  including  the  cohesion  term. 

Tne  program  determines  the  factor  of  safety  of  the  wedge  by  an  incremental 
application  of  the  loads  and  predicts  the  motion  of  the  wedge  in  time.  Rotational 
effects  are  fully  taken  into  account. 
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Fioure  2.  I"  snows  tne  geometry  cf  a  wedge  considered  ir.  this  example. 

Two  cases  were  examined.  In  the  first,  ioints  exhibit  only  strength  that  is 
proportional  to  normal  stress,  i.e.,  there  is  friction  but  no  cohesion.  In  tne 

second,  toints  exhibit  strength  tnat  is  inoepenaent  of  normal  stress,  i.e.,  tnere 

is  conesion  but  no  friction.  Figure  2.16  shows  the  loads  acting  on  tne  Dase  of 
She  wedae  and  Figure  2.19  tne  elements  into  which  the  base  of  the  weage  is 
divided. 

Tne  results  from  the  first  case  (friction  only)  are  presented  ir.  Figure  2.20 
as  a  plot  of  displacements  of  the  wedge  versus  time.  There  is  no  rotation,  as  is 
always  the  case  if  the  external  loads  acting  on  the  wedge  can  be  reduced  to  a 

single  force  and  if  there  is  no  cohesion.  This  is  so  because  shear  stress,  shear 

def ormability  and  shear  resistance,  being  proportional  to  normal  stress,  vary  in 
the  same  way  over  the  failure  surface,  thus  resulting  in  a  translational  motion  of 
the  wedge. 

The  results  from  the  second  case  (cohesion  only'  are  shown  in  Figure  2.21 
as  plots  of  the  displacements  cf  the  wedge  versus  time.  In  this  case 
tnere  is  rotation  of  the  wedge,  because  shear  resistance  is  uniform  over  the 
failure  surface  wnereas  shear  stresses  are  net.  Note,  however,  that  under 
different  external  loading  conditions,  rotation  effects  may  be  present  even  in  tne 
absence  of  cohesion. 

In  Doth  cases  the  loads  were  applied  instantaneously  at  t  =  C  .  Under  a 
very  slow  incremental  application  of  the  loads  (equivalent  to  ignoring  time 
effects),  both  cases  have  factors  of  safety  equal  to  C.47. 

The  wedge/biock  analysis  developed  in  this  research  makes  it  possible  tc 
represent  the  Dehavior  cf  geints,  the  displacements  cf  rock  wedges  and  the 
evolution  of  displacements  in  time.  Kith  these  features  it  is  possible  tc  predict 
and  assess  stability  as  well  as  tne  state  cf  incipient  failure. 
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Figure  3.20  Displacement  Time  Plots  (Friction-only  case) 

There  is  no  rotation.  All  elements  simultaneously 

REACH  6j  AT  T  =  0.88T  AND  fn  AT  T  =  A.18T  (AT 
THAT  TIME  FAILURE  OCCURS), 


21  Displacement  Time  Plots.  The  larger  linear 

PART  OF  THE  PLOTS  IS  DUE  TO  THE  FACT  THAT 
WHEN  TAN  <P  j  =  TAN<Pn  (SEE  EQUATION  3.14). 

Failure  occurs  at  t  =  16. 8T. 


4 .  Fracture  Propagation  in  Jointed  Media 
4 . 1  Background 

Displacements  along  or  normal  to  joint  surfaces  are  major  contributors  to 
rock  mass  de forma bility.  Similarly,  and  as  has  been  shewn  above,  many  failure 
modes  of  rock  masses  involve  displacement  along  joints.  Nevertheless,  rock  mass 
deformability  and  failure  (strength)  can  also  be  caused  by  propagation  of  new 
fractures.  Such  fractures  can  propagate  between  existing  joints  located  in  the 
same  plane  or  they  can  propagate  from  a  joint  in  one  plane  to  one  in  another 
plane.  Various  studies  of  this  behavior  have  been  conducted.  For  instance 
experimental  studies  on  model  rock  masses  show  the  different  types  of  failure 
including  propagation  of  new  fractures  (e.g.  Einstein  and  Hirschfeld  (1973). 

Lajtai  (1969)  shewed  that  shearing  under  low  normal  stresses  usually  leads  to 
tensile  fracturing.  This  approach  was  incorporated  in  stochastic  rock  mass  failure 
models  of  the  MIT  rock  mechanics  group  (see  e.g.  Einstein,  et  al.  ,  1981)). 

The  intent  of  the  research  reported  here  was  to  extend  the  preceding  work. 
Specifically  concepts  and  approaches  of  fracture  mechanics  were  used  to  model  the 
propagation  of  new  fractures  from  existing  joints*.  Also,  appropriate  numerical 
methods  were  developed  for  modeling  existing  joints,  the  initiation  and 
propagation  of  new  fractures, and  the  remainder  of  the  rock  mass  and  its 
boundaries.  These  methods  can  be  combined  with  the  stochastic  rock  joint  system 
models  to  represent  fracturing  of  rock  masses;  initial  work  on  this  last  aspect 
was  also  completed. 

*  Note  that  in  the  following  text  the  term  fracture  is  used  most  of  the  time 
rather  than  joint.  This  is  done  to  indicate  the  generality  of  the  procedure  which 
is  applicable  to  any  "fracture  like"  discontinuity  in  rock. 
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. 2  Hybridized  Displacement  Discontinuity  -  Indirect  Boundary  Element  Method 
The  Displacement  Discontinuity  Method  uses  only  Displacement  Discontinuity 
Elements  and  the  Indirect  Boundary  Element  Method  uses  Stress  Discontinuity 
Elements.  Bv  incorporating  the  two  types  of  elements  in  the  same  numerical 
procedure,  the  two  methods  are  effectively  hybridized.  The  Indirect  Boundary 
Element  Method  can  model  a  continuum  better  than  the  Displacement  Discontinuity 
Method  but  the  latter  method  can  model  fractures  better.  By  modelling  the  external 
boundaries  with  Stress  Discontinuity  Elements,  the  advantages  of  both  methods  are 
utilized. 

It  has  to  be  mentioned  that  both  methods  and  thus  the  hybridized  method  are 
based  on  assumptions  of  linear  elasticity.  For  the  intact  rock  this  is  a 
simplified,  but  usually  justified,  assumption.  For  fracture  propagation  the 
assumption  of  Linear  Elastic  Fracture  Mechanics  (LEFM)  can  be  used  as  a  first 
approximation,  but  should  often  be  replaced  by  more  appropriate  approaches  (see 
Section  4.2.2) . 

4.2.1  Details  on  Elements 

Analytical  influence  functions  (induced  stresses  and  displacements)  of  the 
Displacement  Discontinuity  Elements  and  the  Stress  Discontinuity  Elements  have 
been  aerived. 

A  Displacement  Discontinuity  Element  has  a  positive  and  a  negative  surface 
which  are  separated  by  an  arbitrarily  small  distance  in  the  undeformed  state  (Fig. 
4.1a).  A  separation  and/or  a  slip  can  be  applied  to  the  two  surfaces  to  model 
(part  of)  a  fracture  (Fig.  4.1b,c).  The  element  shown  in  Fig.  4.1  has  a  constant 
displacement  discontinuity  along  its  axis  and  is  called  a  Constant  Displacement 
Discontinuity  Element.  Elements  with  sgue.re  root,  linear  and  parabolic 
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distributions  of  displacement  discontinuities  are  considered.  Each  of  tneir 
influence  functions  are  derived  analytically  by  integrating  tne  induced  stresses/ 
displacements  due  to  the  infinitesimal  dislocations  along  the  element  axis. 

Except  for  the  parabolic  elements,  these  influence  functions  have  been  checked 
with  published  results  and  results  by  numerical  integration. 

A  Stress  Discontinuity  Element  has  a  positive  and  a  negative  surface  which 
are  separated  by  an  arbitrarily  small  distance  in  the  undeformed  state  (Fig. 4. 2a). 
Normal  and/or  shear  stresses  can  be  applied  to  the  two  surfaces  to  model  external 
loads  ( Fig . 4 . 2b , c ) .  As  a  result  of  the  applied  stresses,  ^umps  in  the  induced 
stresses  across  the  element  axis  occur  and  are  termed  stress  discontinuities. 
Elements  with  a  constant  and  a  linear  distribution  of  stress  discontinuities  are 
considered.  Tne  influence  functions  of  the  Linear  Stress  Discontinuity  Element 
have  beer,  derived  analytically  using  Kelvin's  point  force  solution  for  plane 
strain. 

4.2.2  Fracture  Propagation  Criteria 

For  a  closed  fracture  unoer  compression  and  shear,  uneven  shear  stresses 
exist  at  the  fracture  surfaces.  Tne  distribution  of  the  shear  stresses  depends 
both  on  tne  distriDutions  of  tne  normal  stresses  and  the  slip  across  tne  fracture 
surfaces.  Consequently,  it  is  douDtful  wnetner  LEFM  assumptions  (i.e.,  a  slip 
distribution  proportional  to  tne  square  root  distance  from  the  tip)  are  valid  near 
the  fracture  tip.  Tnus  a  more  realistic  element  type  (the  Parabolic  Displacement 
Discontinuity  Element)  which  ooes  not  give  infinite  stresses  at  the  tip  is  used. 
The  stress  state  at  tne  tip  can  De  calculated  and  tne  Maximum  Tensile  Stress 
criterion  is  used  for  modelling  fracture  propacatior. ,  i.e.  tne  fracture 
propagates  ir,  that  radial  direction  where  tne  hoop  stress  reacnes  tne  tensile 
strength  of  tne  material. 
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(a)  undeformed  state 
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(d)  applied  normal  stress 


(c)  applied  shear  stress 


Figure  4,2  Stress  Discontinuity  Elemen 
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For  an  opened  fracture  it  is  assumed  that  no  stresses  act  on  tne  fracture 
surfaces.  LEFM  conditions  are  assumed  near  tne  fracture  tip  and  tne  stress 
intensity  factors  are  calculated.  Tne  Maximum  Tensile  Stress  Factor  criterion  is 
used  for  modeling  propagation. 

4.3  Implementation 

A  Boundary  Element  program  FROCK  (Fractured  Rock)  utilizing  the  two  types  of 
elements  discussed  in  4.2  has  been  developed.  It  has  the  following  advanced 
features : 

(1)  Any  common  type  of  linear  elastic  boundary  conditions  can  oe  modelled. 
These  include  stress,  displacement,  mixed  (e.g.  roller)  and  spring 
boundaries . 

( 2 )  Automatic  treatment  of  non-linear  responses  due  to  non-linear 
stress-slip  relations  (e.g.  slip-weakening)  of  closed  fractures, 
opening/closing  of  factures,  and  fracture  propagation. 

(3)  Since  analytical  influence  functions  of  the  elements  are  used,  only 
single  precision  is  sufficient  for  the  Fortran  program  FROCK. 

(4)  Optional  acceleration  of  non-linear  analyses. 

(5)  Optional  iterations  to  achieve  better  accuracy.  Usually  a  single 
additional  iteration  is  sufficient. 

(6)  Modulized  input  commands  -  the  user  can  modify  different  parts  of  the 
problem  at  will  during  the  analysis. 

(7)  Mesh  generation  aids  to  facilitate  input. 

(8)  Tne  problem  status  can  be  stored  and  further  analyses  can  De  done  at  a 
later  time. 

(°  Realistic  modelling  of  ir,  situ  conditions  is  obtained  oy  incorporating 
Gravity  stresses  as  initial  stresses. 


4.4  Theorems 


It  was  found  tnat  the  influence  functions  of  the  Displacement  Discontinuity 
Elements  do  have  a  numoer  of  interesting  characteristics.  The  induced  stresses/ 
displacements  due  to  the  opening  up  of  an  element  are  approximately  proportional 
to  the  total  opening  or  "expansion"  at  distances  far  away  from  the  element  and  are 
independent  of  the  opening  shape.  Thus  the  "Theorem  of  Equivalent  Expansions"  was 
formulated  which  may  have  important  potential  applications  such  as  the  estimation 
of  the  expansion  of  underground  fractures  by  surface  or  near-surface  measurements. 
Similar  observations  were  made  on  the  relative  slip  of  the  elements.  For  Stress 
Discontinuity  Elements  the  induced  stresses/aisplacements  at  large  distances  are 
found  to  be  approximately  proportional  to  the  total  applied  force  on  the  element 
and  to  be  independent  of  the  pressure  distribution  on  the  element.  The  "Theorem 
of  Equivalent  Forces"  was  formulated  which  is  in  fact  similar  to  Saint  Venant 1 s 
principle . 

4.5  Case  Studies 

Many  case  studies  were  conducted  to  demonstrate  the  versatility  of  the 
program.  FROCK  and  to  check  the  validity  of  the  results.  These  cases  are: 

(1)  Opened  fractures  in  the  infinite  medium  under  tension  and/or  shear 
stresses.  Various  configurations  have  been  tried  (Fig. 4.3). 

The  stress  intensity  factors  calculated  by  FROCK  are  compared  with 
published  results  whenever  available. 

(2)  Internally  fractured  rectangular  plate  under  tension  (Fig.  4.4). 
Different  element  types  and  boundary  conditions  are  used  in  a  number  of 
cases.  The  K- ' s  calculated  are  usually  within  15%  of  published  results. 


(a)  1  fracture 


(b)  stepping  fractures 


c)  parallel  fractures 


(d)  unequal  fractures 


colinear  fracture 


(f)  2 


inclined  fractures 


Figure  4,3  Fracture  Configurations  Considered 


Direct  Shear  Test.  A  Direct  Shear  Test  on  a  rock  joint  which  has  a 


non-linear  relation  between  the  shear  stress  t  and  the  slip 
displacement  6  is  simulated  (Fig.  4.5).  The  diplacement  jump  across 
the  joint  is  shcvm  by  plotting  the  horizontal  displacement  ux  along  a 
vertical  section  (Fig.  4.6).  The  shear  stresses  along  the  joint  may  not 
be  uniform  and  are  highest  in  the  center  part  of  the  joint.  Fig.  4.7  is 
a  plot  of  the  shear  stress  along  the  joint  when  the  applied  horizontal 
displacement  is  6  x  10-4  (consistent  units)  at  the  top. 

Fracture  propagation  in  a  plate  under  compression:  The  propagation  of  a 
fracture  in  a  rectangular  plate  under  compression  is  simulated  (Fig. 
4.8).  The  propagation  trajectories  from  the  initially  straight  fracture 
seem  to  be  consistent  with  experimental  observations  and  theoretical 
s  tipula  tions. 

Fracture- fracture  interaction  and  propagation  under  vertical  tension  in 
an  infinite  medium:  Different  configurations  are  considered  (Fig.  4.9). 
For  stepping  fractures  with  small  or  no  overlap  the  general  pattern  is 
that  initially  the  inner  tips  have  higher  stress  intensities  than  the 
outer  ones  and  propagate  first.  After  a  while  the  propagations  may  or 
may  not  cause  intersection  of  the  two  fractures.  Then  the  outer  tips 
start  to  propagate  and  approach  the  horizontal  direction.  Further  on 
the  outer  tips  open  up  essentially  in  Mode  I  with  Kj's  approximately 
equal  to  those  of  a  horizontal  crack  with  width  equal  to  the  horizontal 
distance  beween  the  two  fractures  tips  (Fig.  4.10). 

Fracture  propagation  in  a  plate  under  vertical  tension:  The 

propaga tion  of  an  inclined  fracture  under  vertical  tension  is  simulated 


(Fig.  4.11).  The  trajectories  compare  well  with  published  results  from 


equivalent  horizontal  fracture 


Figure  4,10  Fracture-Fracture  Interaction  and 
Propagation 
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a  fatigue  experiment  on  a  titanium  alloy  (Fig.  4.12).  The  propagations 
of  2  stepping  fractures  in  a  plate  is  also  simulated  with  satifactory 
results . 

7)  Behavior  of  a  joint  in  infinite  medium  under  compression:  Details  of 
the  behavior  of  an  inclined  joint  in  an  infinite  medium  under  vertical 
compression  were  examined,  particularly  the  developements  of  sliding 
and  tension  at  the  tips  with  increasing  compressive  stress. 

8)  Fracture  propagation  in  a  plate  under  compression:  The  effects  of 
gravity,  horizontal  and  vertical  compressive  stresses,  and  free  ground 
surface  on  propagation  were  examined. 

9)  Stability  of  an  excavated  slope:  The  possible  propagation  of  new 
fractures  from  joints  in  an  excavated  slope  leading  to  slope  failure  was 
examined.  At  first  a  slope  without  joints  was  considered  and  its  stress 
distributions  examined,  followed  by  a  slope  with  deterministic  joints. 
Finally  the  uncertainty  of  the  joint  geometry  was  considered  to  some 


extent. 


-  experimental  trajectory 

•  .  •  •  •  output  from  FROCK 

Figure  4,12  Comparison  of  Numerical  and  Experimental  Results 
(Tension  Test  on  Titanium  Alloy) 
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5.  Conclusions  and  Furtner  Research 


The  research  reported  here  has  produced  a  number  of  significant  contributions 
in  the  area  of  rock  mass  behavior: 

•  Joint  geometry  models.  Two  and,  particularly,  three-dimensional 
stochastic  joint  system  models  and  associated  computer  codes  have  been 
created.  They  represent  many  actual  rock  mass  geometries,  relying  on 
standard  statistical  sampling  methods,  and  providing  the  basis  for  flow 
and  fracture  analysis. 

o  Flow  through  pointed  media.  The  stochastic  joint  system  models  have 
been  combined  with  pipe  network  and  Finite  Element  flow  models  to 
represent  flow  through  jointed  media  in  two  and  three  dimensions.  The 
associated  computer  codes  have  been  developed.  Realistic  flow 
analysis  in  rock  is  thus  possible  which  is  needed  for  engineering 
design  (cut  slopes,  tunnels,  dams,  waste  storage)  and  for  hydrogeologic 
applications.  Use  in  resource  extraction  can  be  foreseen. 

•  Simplified  expressions  for  joint  geometry  and  hydraulic 
characteristics  of  jointed  rock  masses  have  been  created.  Initial 
parametric  studies  show  that  "joint  area  per  volume"  (persistence 
measure  P32)  and  "connectivity"  are  simplified  measures  which  can 
represent  rock  mass  behavior  and  thus  replace  the  complete  models. 

•  Parametric  studies  with  the  flow  in  jointed  media  models  have 

shown  that  joint  system  size  relative  to  the  domain  (scale  effect)  and 
anisotropy  (isotropy)  of  the  joint  systems  have  significant  effects. 

Also  three  dimensional  approaches  lead  to  results  that  differ 
significantly  from  two  dimensional  approaches.  Finally,  joint  aperture 
does  not  seem  to  have  always  the  overriding  effect  that  was  stipulated 
previously . 
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•  By  combining  the  distinct  element  method  and  an  analysis  method  for  flow 
in  joints  of  variable  aperture  a  f low/def ormation  analysis  for  jointed 
rock  masses  was  created.  Associated  interactive  computer  programs  have 
been  developed.  The  method  couples  deformation  and  flow,  i.e.  the  effect 
of  seepage  pressure  changes  joint  apertures,  and  changes  in  joint 
apertures  affect  flow.  This  analysis  method  is  an  essential  tool  for 
designing  slopes,  tunnels,  wells  and  other  structures  in  jointed  rock.  It 
is  also  essential  for  interpreting  monitoring  observations  on  water 
pressure  and  movements  in  these  engineering  projects, 
o  Extensive  parametric  studies  with  the  flow/deformation  analysis  have  been 
conducted  for  slopes  in  jointed  rock,  injection/withdrawal  wells  in 
jointed  rock,  and  dams  on  jointed  rock  (with  and  without  grout  and 
drainage  curtains ) .  The  results  show  that  coupling  is  significant  if  the 
joint  apertures  vary  within  the  rock  mass.  Ignoring  the  effect  of 
coupling  may  lead  to  unconservative  designs, 
o  a  wedge/block  analysis  and  computer  code  was  developed  with  which  the  com¬ 
plete  movement  of  wedges  (blocks)  in  translation  and  rotation  can  be 
modeled,  and  which  can  represent  the  prefaiiure  time  dependent  behavior. 

®  A  numerical  method  for  fracturing  of  jointed  media  has  been  developed 
based  on  displacement  discontinuity  elements  of  different  types  and 
boundary  elements.  Bounded  and  unbounded  domains  can  be  represented,  as  can 
fracture  in  mode  I  and  II,  and  sliding  of  joints  or  newly  propagated 
fractures  in  a  slip  weakening  mode.  This  analysis  method  is  essential  to 
represent  the  deformation  and  possible  failure  of  rock  masses  in  which  a 
combination  of  joint  sliding  and  fracture  of  intact  rock  takes  place.  This 
method  has  been  combined  with  the  joint  geometry  models. 
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•  The  fracture  model  for  pointed  media  produces  tne  experimentally 
observed  propagation  of  new  fractures  into  a  direction  parallel  to  tne 
major  principal  stress.  It  also  produces  fracturing  in  Detween  joints 
and  between  joints  and  free  surfaces  which  are  important  in  rock  mass 
failure.  Of  particular  interest  are: 

Opening  of  fractures  (joints)  in  tension  for  both  infinite  and 
bounded  fields  produces  stress  intensity  factors  which  are  close 
to  the  theoretically  predicted  ones. 

Fractures  (joints)  in  compressed  stress  fields  propagate  into  the 
direction  of  the  major  principal  stress  as  theoretically  predicted 
and  experimentally  shown.  Interaction  of  several  fractures 
(joints)  and  different  inclinations  to  each  other  and  to  the 
stress  field  can  be  modeled.  The  sequence  of  propagation  is 
intricate  and  switches  from  one  end  of  the  existing  fracture 
(joint)  to  the  other. 

-  A  number  of  tneorems  expressing  the  effect  of  existing  joints 
(fractures)  on  the  behavior  of  the  jointed  rock  mass  under 
different  stress  fields  have  been  developed. 

Although  the  research  results  are  significant  and  represen\  major 
advances  further  work  is  necessary  in  a  number  of  areas: 

•  The  joint  geometry  models  need  to  be  thoroughly  validated  in  the  field 

•  Modelling  of  nonsaturatec  flow,  multiphase  flow  and  transport  would 
be  of  interest. 

o  The  simplified  measures  need  to  be  verified  for  flow  application?  a:., 
further  developed  for  deformation  (failure)  applications. 
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•  Extension  of  the  coupled  flow  deformation  analysis  to  three  dimensions 
and  other  block  shapes  is  desirable. 

•  Implementation  of  the  combined  stochastic  geometry  and  fracture  models 
with  field  or  laboratory  verification  is  needed. 

•  The  use  of  non-linear  fracture  mechanics  needs  to  be  considered  to 
correctly  model  fractures  approaching  free  surfaces  or  other  joints. 
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7 .  Publications 

Publication  of  research  results  is  accomplished  throuqh  this  final  report 
series  and  a  number  of  papers  in  journals  and  symposia/conferences. 

Published  Papers: 

"Coupled  Deformation/Flow  Analysis  with  the  Distinct  Element 
Method"  by  J.  Kafritsas,  M.  Gencer.  H.  H.  Einstein,  Proc . 

Symp.  25th  U.S.  Symp.  on  Rock  Mechanics,  1984. 

"Methods  of  Artificial  Intelligence  in  Rock  Mechanics"  by 
K.  D.  Dershowitz,  H.  H.  Einstein,  Proc.  Symp.  25th  U.S.  Symp. 
on  Rock  Mechanics,  1984. 
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"A  Three-Dimensional  Model  for  Flow  of  Fractured  Rock"  by  W.  S. 
Dershowitz,  B.  M.  Gordon,  J.  Kafritsas,  H.  Herda,  Proc .  17th  Int. 
Congress  of  the  Int.  Assoc,  of  Hydrology,  1985. 


Papers  submitted  for  Publication: 

«  "Conceptual  Models  for  Characterization  of  Rock  Joint  Geometry" 
by  W.  S.  Dershowitz,  H.  H.  Einstein,  submitted  to  Rock  Mechanics 
and  Rock  Engineering. 

o  "Scale  Effects  in  the  Hydrology  of  Jointed  Rock  Masses"  by  W.  S. 
Dershowitz,  H.  H.  Einstein,  submitted  to  Water  Resources 
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Planned  Papers : 

•  Paper  on  parametric  studies  with  2-D  and  3-D  joint  flow  models. 

•  Paper  on  simplified  measures  for  expressing  joint  geometry  and 
flow  on  jointed  media. 

•  Paper  on  parametric  studies  with  flow/deformation  analysis. 

•  Paper  on  wedge /block  analysis  with  rotation  and  time  dependent 
behavior . 

•  Two  papers  on  fracturing  of  jointed  media,  one  on  fracture 
mechanics  principles  and  numerical  approaches,  one  on  results. 
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